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Fundamental parameters of RR Lyrae stars from 
multicolour photometry and Kurucz atmospheric models 
II. Adaptation to double- mode stars 
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1 INTRODUCTION 



As described in the first paper of this series (|Barczall2O10l . 
hereafter Paper I) a new method can be used to deter- 
mine fundamental parameters of RR Lyrae (RR) stars using 
broad-band optical photometry and the conservation laws of 
mass and momentum in the pulsating atmosphere. 

The first version of the method (|Barczal 12001 |2006| 1 
used the law of momentum conservation in the frame of a 
uniform atmosphere approximation (UAA), that is, the pul- 
sation of the atmosphere is taken into account as if the at- 
mosphere were a rigid shell. The available Johnson- Cousins 
UBV(RI)c photometries of SU Dra and T Sex (|Barczal 
l2002i , l2006t ) were processed as examples because these uni- 
formly cover the whole cycle of pulsation and allow a solu- 
tion of the Euler equation of hydrodynamics for the mass 
M a of the star and distance d to it. (The subscript 'a' in- 
dicates that this mass is a dynamical mass derived from an 
analysis of the motion of the atmosphere.) 

In Paper I, an extended hydrodynamic treatment, in 
which the UAA is dropped, was reported. The following two 
main steps were involved. 

(i) Assuming a perfect spherical symmetry of the pulsa- 
tion, the photometric quantities (colour indices and bright- 
ness) were converted to time-dependent physical quantities 
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ABSTRACT 

Our photometric-hydrodynamic method is generalized to determine fundamental pa- 
rameters of multiperiodic radially pulsating stars. We report 302 U BV(RI)c Johnson- 
Kron-Cousins observations of GSC 4868-0831. Using these and published photometric 
data of V372 Ser, we determine the metallicity, reddening, distance, mass, radius, equi- 
librium luminosity and effective temperature. The results underline the necessity of 
using multicolour photometry, including an ultraviolet band, to classify the subgroups 
of RR Lyrae stars properly. Our U observations might reveal that GSC 4868-0831 is 
a subgiant star pulsating in two radial modes and that V372 Ser is a giant star with 
size and mass of an RRd star. 
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stars: individual: GSC 4868-0831 - stars: individual: V372 Ser - stars: variables: RR 
Lyrae. 



(effective temperature T e , effective gravity g e and angular 
radius 1?) using the computed colours and fluxes of the AT- 
LAS models of lKuruczl |l997l ). 

(ii) The physical quantities were introduced in the differen- 
tial equations expressing the laws of mass and momentum 
conservation during the pulsation. Particular solutions were 
given to describe the motion of the pulsating atmosphere 
in the gravity field of the star. The two time-independent 
parameters of the solutions - the mass A1 a of the star and 
the distance d to it - were determined. 



Because the ATLAS models apply to the atmosphere of 
non-variable stars, quantitative photometric and hydrody- 
namic conditions (Conditions I and II in Paper I, hereafter 
C (I) and C (II) , respectively) were formulated for the applica- 
bility of the quasi-static atmosphere approximation (QSAA) 
in order to find the time intervals of the pulsation when dy- 
namical phenomena have a negligible effect on the colours 
and brightness, (i.e. the structure and colours of the atmo- 
sphere are identical to those of a selected ATLAS model). 

A summary of the conditions is as follows. C (I) is satis- 
fied if the difference of the continuum fluxes of the observed 
and selected ATLAS model does not exceed the error of the 
observation in the optical spectrum covered by the colours 
U,...,Ic- G^ 11 ' is satisfied if the acceleration in the atmo- 
sphere is equal to the inst antaneous 'effective gravity' g e (t) 
l|Ledoux fc Whitnevl[l960l ) of the selected ATLAS model. 

The UBV(RI) C photometry of the RRab star SU Dra 
was used to demonstrate that the extended method is a 
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viable alternative to determine the fundamental parameters 
of RR stars. The atmospheric metallicity [M] , the reddening 
E(B — V) towards the star, d and .Ma were determined from 
phases when the conditions of the QSAA were satisfied. 

Double-mode (DM) RR (RRd) stars pulsate in two 
radial modes simultaneously. Their importance for stellar 
pulsation theory is obvious because DM pulsation offers 
a unique possibility to determine fundamental parameters 
such as mass Ai p and luminosity L p from the Petersen dia- 
gram, (i.e. from frequencies that are accessible by observing 
the brightness variation over a sufficiently long time scale). 
M p and L p can be compared with the mass A^ev and lu- 
minosity L cv derived from stellar evolution theory. The as- 
sumptions M ev = M P , L DV = L p plus some colour informa- 
tion have gi ven, for example, t he fundamental parameters 
of BS Com jDekanv et allhoosl ). 

In this paper, we use our combined photometric- 
hydrodynamic method to determine the fundamental pa- 
rameters d and Ala, the approximate position in a theoret- 
ical Hertzsprung- Russell diagram (HRD), the radius varia- 
tion, the reddening and the metallicity of the DM pulsators, 
GSC 4868-0831 and V372 Ser. We also describe the kine- 
matic behaviour of the pulsating atmosphere in our lim- 
ited hydrodynamic treatment. The accuracy of the light and 
colour curves cannot be enhanced by folding, and therefore 
we give a refinement of the technique described in Paper I. 
The method can be applied to RR stars with any number of 
periods (one or > 2) if there are photometric observations 
available in sufficient numbers. However, an application to, 
for example, multiperiodic 5 Set stars with small amplitudes 
would allow the determination of the A^arf -2 only because 
the hydrodynamic status of the atmosphere has only a small, 
non-radial variation. 

A4 a derived here provides a mass value from a com- 
pletely different astrophysical input in comparison with M ev 
or Alp. Consequently, A4 a can be an independent check for 
evolution and pulsation theory of RR stars. To the best of 
our knowledge, we are the first to attempt to determine the 
distance and mass, etc. of DM pulsators with an astrophys- 
ical method using only the motion of the atmosphere. 

In comparison with the Baade-Wesselink (BW) method, 
the main advantage of our method is that more output is 
obtained for less observational input. Spectroscopic observa- 
tions are not necessary at all, and consequently our method 
can easily be applied to faint stars. A BW solution has not 
been found for RRd stars in the literature. Perhaps this 
can be explained by the problems arising fro m the faint- 
ness (V > 10.5 mag for the known DM stars, IWilsl 120061 . 
ISzczygiel fc Fabrvczkvl 120071 ) and multiperiodic character 
of DM pulsation: simultaneous observation of light, colour 
curves and spectroscopy of faint stars would be necessary 
over days. Furthermore, problems are encountered with the 
precise determination of the centre-of-mass velocity, a sub- 
stantial point of the BW analysis (Paper I). Our extended 
photometric-hydrodynamic method promises to deliver pa- 
rameters in addition to distance, the main fundamental pa- 
rameter that can be acquired by the BW analysis. 

A recent challenge to pulsation theory originates from 
the MOST satellite, which found frequencies of the RRd 
star, AQ Leo, with amp litudes down to the mmag level 
jGruberbauer et al.ll2007l ). Our method is an extension of 
the research methodology of DM stars beyond theoretical 



Table 1. Log of the observations of GSC 4868-0831. 



HJD-2 400 000 


No. of frames 


Telescope 


54822. 5555-. 5868 


35 


RCC 


54829.4534-.6300 


225 


RCC 


54830.4442-.4950 


75 


RCC 


54831.4606-.6221 


215 


RCC 


54832. 4418-.6122 


205 


RCC 


54863.4483-. 7015 


265 


IAC80 


54871. 4757-.6242* 


70 


IAC80 


54873.3661-.5941 


235 


IAC80 


54874.4437-.4620* 


20 


IAC80 


54876.4205-.6108* 


165 


IAC80 


* Epoch of the tie-in 


observations 





and empirical methods using only frequencies and ampli- 
tudes, for which data can be obtained from a single-band 
time series. 

The observations, standard UBV(RI)c magnitudes of 
the variables and some field stars are reported in Section 2. 
The metallicity and reddening of the variables and compar- 
ison stars are given in Section 3. As a by-product, T a ,logg, 
angular radius $ are also determined for the comparison 
stars. In Section 4, we describe some technical details beyond 
those reported in Paper I. The results are presented here for 
the brightest DM pulsators GSC 4868-0831 and V372 Ser, 
and an insight is given into the kinematics of their atmo- 
spheres. We give a discussion and our conclusions in Sections 
5 and 6, respectively. In Appendix we describe the publicly 
available program package BBK0 which can be used to de- 
termine the fundamental parameters from the photometric 
input. 



2 THE OBSERVATIONS AND REDUCTION 

The two brightest DM pulsators are the subject of the this 
paper. The observational data of GSC 4868-0831 were col- 
lected with the IAC8C0 telescope of the Teide Observa- 
tory and the 1-m RCC telescope mounted at Piszkesteto 
Mountain Station of the Konkoly Observatory. The tech- 
nical description of the CCD cameras and the details 
of the observations and r eductions are identical to those 
of iBenko fc Barczal (|2009T ) . We used the standard lRAi{f| 
tasks for the reductions. Transformation into the standard 
U BV(RI)c syste m was done by using the equatorial stars 
of lLandoltl d 198St). The observations of V372 Ser have been 
published in Ben ko fc Barczal (|2009h . A log of the observa- 
tions of GSC 4868-0831 is given in TableQ] typical exposure 



The program package is available from 
http://www.konkoly.hu/staff/barcza/pub.html It is composed 
of tables extracted from http://kurucz.harvard.edu/grids.html, 
FORTRAN source codes and a manual. 

2 The 0.82m IAC80 Telescope is operated on the island Tenerife 
by the Instituto de Astrofisica de Canarias in the Spanish Obser- 
vatorio del Teide. 

3 I RAF is distributed by the National Optical Astronomical Ob- 
servatory, operated by the Association of Universities for Re- 
search in Astronomy Inc., under contract with the National Sci- 
ence Foundation. 
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GSC 4868-0831 



V372 Ser 



£ 10.2 

> 



1 %\ 



Phase 

Figure 1. V light curves from our observations. The upper and 
lower panels are folded with Pi and Pq, respectively. 



Table 2. Result of the photometry giving the average values 
of the two variables and some comparison stars in the field of 
GSC 4868-0831. The magnitude errors of the comparison stars 
are 0.008, 0.011, 0.007, 0.009, 0.009 mag in V ,U - B.B - V,V - 
Rc,V — Ic, respectively. 



ID 


V 


U -B 


B - V 


V -Rc 


V -I C 


V372 Sort 


11.350 


0.000 


0.380 


0.256 


0.524 


GSC 4868 












-0831* 


10.769 


-0.087 


0.357 


0.215 


0.478 


-0063 


12.504 


-0.12 


0.468 


0.287 


0.583 


-0436 


12.637 


0.220 


0.648 


0.365 


0.705 


-0779* 


12.534 


0.092 


0.554 


0.331 


0.641 


-0860* 


12.579 


0.045 


0.496 


0.278 


0.519 


-1089* 


12.862 


0.073 


0.594 


0.353 


0.673 



Notes. 

* ■ The magnitude averaged value from TV = 529 observations. 
The magnitude averaged value from 280 observations. 
The stars are visible only in the frames taken with IAC80. 



times were 180, 40, 20, 6, 10 s for U, B, V, (R, I)c, respec- 
tively. The folded V light curves are plotted in Fig. [T] 

The result of the tie-in observations for the compar- 
ison and check s tars o f V372 Ser is given in table 2 of 
iBenko fc Barczal (120091 1 and those for GSC 4868-0831 are 
given in Table [2] Additionally, we give the magnitude av- 
eraged V and colour indices of the variables obtained from 
our observations. They characterize the variables because 
the distribution of the observations is quasi-random over a 
long enough time. 



2.1 GSC 4868-0831 

The X — U,B,V, Rc and Ic magnitudes were obtained 
by differential photometry from the instrumental Ax = 
Au,A6,... magnitude differences of GSC 4868-0831, GSC 
4868-0063 and GSC 4868-0436 as follows. The coefficients 
cx ~ 1 and zero points Cq « of the transformation equa- 
tions 



^GSC4 



= ^comp + cx Ax + c 



(1) 



were determined from all frames separately for the telescopes 
IAC80 and RCC, respectively. An average was computed 
from the comparison stars GSC 4868-0063, GSC 4868-0436. 
Next, the magnitudes U, B, Rc and Ic were interpolated to 
the epoch of V and the colour indices were computed. These 
are the data that form the basis for our analysis 



2.2 V372 Ser 



/ X) 

We carefully revised the zero points A nn of the light 
curves of V372 Ser given in table 4 of iBenko fc Barczal 
(2009) because it is particularly important to have 
the colour indices in the standard system. An error in 
the zero point of the magnitude scales was removed 
which resulted in the shifts AX of the amplitudes 
Ann' 1 ^ = A ( o ] + AX in Eq. (2) of IBenko fc Barczal 
(2009) where AU = +0.020, AB = +0.010, AV = +0.049, 
AR C = +0.029, Ale = -0.017. The revised colours 
U, B , Rc , Ic were interpolated to the epoch of V and the 
colour indices were obtained by subtractions. The averages 
of the revised colour indices are now in better agreement 
with those of iGarcia-Melendo. Henden fc Gomez-Forreladl 
i|200ll) . The observations of HJD - 2454200 = 
17.3765 20.42 -. 47, 44.34-. 38, 45.3464, 45.3808 
|Benko fc Barczal 2009) were omitted because of 
poor weather conditions. Our final list contains 
V,B — V,U — B,V — Rc, V — la points for N = 529 
epochs. These are available in electronic forrr[f|. 



3 METALLICITY, REDDENING 

The physical quantities, especially log go, derived from pho- 
tometry depend on [M] and E(B — V). [M] is identica l with 
the parameter of the ATLAS models (|KuruczHl997T ). It is 
given on the solar scale (i.e. [M] = dex for the solar com- 
position). Various methods to determ ine [M] and E(B — V) 
are summarized in iLiu fc Janes] (|l990l ). none of which is ap- 
plicable for our stars except for the upper limit of E(B — V), 
which can be found from the maps of iBurstem fc Heilea 
i 19821) and from the diffuse inf rared backgrounc0 (DIRBE, 
ISchlegel, Finkbeiner. fc Davieslll998l ). The precise values are 
of primary importance because we use g e , T e in the Euler 
equation of hydrodynamics. Therefore, we have to use the 



The 



table containing V, B - V, U - B, V — Rc, 
V — Ic observations is available electronically: 
http://www.konkoly.hu/staff/benko/pub.html Flag d indi- 
cates the observations omitted from the analysis because of poor 
sky conditions. 

5 http:/ /www. konkoly.hu/staff/bcnko/pub. html 

6 http : //irsa. ipac . caltech.edu/applications/DUST 
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GSC 4868-0831 



U- 



V372 Ser 




0.5 



0.3 0.4 0.5 



0.6 

V-L 




Figure 2. Colour-colour diagrams (U - B)-(B - V) and (B - V)-(V - Ic)- The panels on the left and right show GSC 4868-0831 and 
V372 Ser, respectively. The solid gree n circles and red crosses denote that CW is satisfied and is not satisfied, respectively. The iso-T e and 
iso-logg e lines of the ATLAS models (Kurucz 1997) are interpolated to E(B — V) and [M] as given in Table|3] The lines with open squares, 
triangles, circles denote logg e = 4.5,3.5,2.5 cms -2 , respectively. The dashed lines from left to right denote T c = 7500, 7000,6500 K. 



photo metric variation method described in lBarcza fc Benko] 
{•21)09.) and Paper I: [M] and E(B - V) are found from the 
best fit of the ob served and th eoretical colour indices of the 
ATLAS models (|Kuruczll 19971) . 

When observing variable stars, we have a sum of the 
random scatter and temporal physical change in the colour 
indices, the effect of the temporal physical change manifests 
itself in increasing A log <? c and AT C . To minimize this source 
of systematic error, TV' 1 ' phases must be selected when the 
photometric condition C' 1 ' of the QSAA is satisfied. The 
averages 

(A log g c ) N m = Alog5e,j, (2) 

i=i 

N m 

(AT C ) N(1) = A^ATe.j (3) 
i=i 

must be minimized as a function of [M] and E(B — V). 

The application is straightforward for the comparison 
stars by setting = 1 and logg c = log<? because their 
atmosphere is static and the observed colour indices have 
a random scatter only. For orientation, E(B — V), [M], 
T e , logg, •& of comparison stars of approximately the same 



colours were determined using our method. The minimal 
value of the scatters Alogg, AT e shows the effect of the 
random error of the colour indices. The averages for the 
different comparison stars suggest that the limits of the ap- 
plicability of the QSAA are AT C < 5 and A log g< 0.04 in 

our observational material. That is, is satisfied in the 
phases where AT e and Alogg e do not exceed these limits. 

The numerical results are (A log 30)39 = 0.022 ± 0.001 
(AT C ) 39 = 8.6 ± 0.5 K for GSC4868-0831, <Alog 5c ) 24 9 = 
0.031 ± 0.002, (AT C >249 = 6.2 ± 0.2 K for V372 Ser. The 
AT C < 5 K limit from the comparison stars is exceeded 
slightly. The excess originates from the inclusion of phases 
into N m , which satisfy A log g c < 0.04 but violate AT C < 5 K 
because of the temporal neighbourhood of an atmospheric 
shock. The results are summarized in Table [3] The atmo- 
spheres of both variables are moderately metal deficient. 

3.1 The observed U - B, B -V, V - I c colour 
indices 

The observed (17 - B)-(B - V), (B - V)-(V - I c ) colour- 
colour diagrams are plotted in Fig. [2] The phases are plotted 
separately when is or is not satisfied (green circles and 
red crosses, respectively). The theoretical colour-colour re- 
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Table 3. Reddening, metallicity, surface gravity, effective temper- 
ature and angular radius of the stars. The errors are AE(B-V) = 
±0.01, and A[M] = ±0.05 for the comparison stars. 



ID 


E[B — V) 


nvii 


1°§ 9e 


T 

± c 


■d x 10 11 




[mag] 


[dex] 


[cms -2 ] 


[K] 


[rad] 


GSC 4868 












-0063 


0.01 


-0.09 


3.98 


6550 


5.05 


-0779 


0.00 


±0.20 


4.07 


6284 


5.31 


-0860 


0.01 


±0.70 


5.20 


6880 


4.92 


-1089 


0.01 


-0.15 


3.08 


6148 


4.90 


-0831+ 


0.008 


-1.05 


3.69* 


6902* 


10.17* 




±0.002 


±0.10 


0.50 


262 


0.32 


GSC 5002 












-0506 


0.005 


-0.10 


4.03 


6247 


8.82 


-0566 


0.01 


-0.05 


3.83 


6228 


4.62 


V372 Ser* 


0.003 


-0.53 


3.24* 


6713* 


8.13* 




±0.003 


±0.05 


0.36 


323 


0.16 



Notes. 

t ATP) = 39. The colour curve segments HJD-2454800=71.6077- 
71.6220,73.5203-73.5919,76.4228-76.5088 were used in Eqs. J2l3t . 
* Averaged values of logg e ,T e ,i? from all observations. In the 
next row, the estimated standard errors of E(B — V), [M] and 
the standard deviations of logg e , T e and $ X 10 11 are given, 
t jVW = 249, the colour curve segments HJD-2454200=17.3765- 
17.5938, 23.3745-23.5858, 42.3402-42.5119, 45.5006-45.6265, 
48.3929-48.4986, 51.4163-51.6057 were used in Eqs. (2\ty . 



lations of the ATLAS models (|Kurucd fl99?f) belonging to 
some characteristic logg c and T c were interpolated for the 
[M] and E(B~ V) values of the variables. These are the lines 
in Fig. [2] Before discussing the hydrodynamic details, we can 
draw some conclusions from the data plotted in Fig. [2] 

• The segregation of the colour indices satisfying C"' is 
clearly seen. However, a remarkable dichotomy is obvious: 
the phases of GSC 4868-0831 satisfying C' 1 ' are concentrated 
in the blue B — V domain, while those of V372 Ser populate 
the red B — V region. 

• Within the observational scatter, the (B — V)-(V — Ic) 
values of V372 Ser are in the domain of 2.5 < logp c < 4.5. A 
considerable number of (B - V)-{V - I c ) pairs of GSC 4868- 
0831, which do not satisfy C^, are significantly below this 
domain indicating logg c S> 4.5. A combination of [M] and 
E{B — V) could not be found that would have shifted every 
pair into a common domain of log<7 e ,T c , [i.e. C™ is more 
strongly violated by the atmosphere of GSC 4868-0831]. 

• As indicated by the difference in log g c — 3.69, 3.24 (Ta- 
ble [3]), the surface gravity of GSC 4868-0831 must be by a 
factor « 3 larger than that of V372 Ser, [i.e. GSC 4868-0831 
is a subgiant rather than a giant star]. 



4 DISTANCE, MASS AND ATMOSPHERIC 
KINEMATICS 

The dynamical equation of the pulsation of an atmosphere 
with spherical symmetry was derived in Paper I from the 
Euler equation of hydrodynamics. Its most convenient form 
is 



where a(r,t) — dv/dt + v(dv/dr), G is the Newtonian grav- 
itation constant. The dynamical correction a^ dyn \r,t) ac- 
counts for the difference between the accelerations of a static 
and a dynamical model atmosphere at time t. The solution of 
the continuity equation for mass conservation in the frames 
of QSAA resulted in a series expansion 



v(r, t) =m - — -^-r ± 
fto ot 



■<R. 



(5) 



This was given in Paper I with detailed explanations for 
the symbols in Eqs. ((4}, ((5}. R is the stellar radius [i.e. 
the radius of approximately zero optical depth], # = R/d, 
the dot denotes a differentiation with respect to t, ho = 
/j,g c (t) /TZT(R, t) is the reciprocal barometric scaleheight at 
r ^ R. 

The velocity profile ((5| was introduced in Eq. ((4]) for 
each epoch j, j = 1, N and 7M a and d were determined in 
four steps. In comparison with the technique in Paper I, the 
refinement of the solution of Eq. Q was rendered possible 
by the large number of the observed points. (N — 280 and 
529 for GSC 4868-0831 and V372 Ser, respectively.) 

(i) The photometric inverse problem [i.e. the conversion of 
the UBV(RI)c observations to log£r e and ho] was solved 
for N observations, as described in Paper I, and the 
phases were found that satisfied C' 1 '. Polynomial fits of de- 
gree 7-10 were calculated for ^}(t),\ogg c (t),ho(R,t). One 
polynomial was sufficient for the whole time interval (of 
length < 0^3) if the atmosphere was in a shock-free state; 
< 3 polynomials were necessary in a strongly shocked phase. 

(ii) Epochs of number n could be selected by differentiation 
of #(t) when v(r,t)/d = constant, and the angu- 
lar acceleration a (ang) (r,t) = a(r,t)/d = i?-i9/i ( 7 1 (a/i /9t)± 
■ ■ ■ ~ over the whole atmosphere because the leading terms 
in v/d and a/d are large in comparison to the rest. We can 
well assume that C^ 11 ' was satisfied in these shock-free in- 
tervals if the atmosphere was in free fall [i.e. a < - dyn '(r, t) ~ 
holds in their vicinity] . This is the descending branch in the 
light curve. Eq. @ reduces, in these epochs, to 

(M a d- 2 ), = G^ 1 g c (U)-& 2 (U) , i = 1, ...,n, (6) 

n is given in Table U These n epochs form a subset of the 
N (1I) epochs satisfying C (II) . 

We remark that a (dyn) (r, t) « might not be expected 
in the minimum and ascending branch of the light curve 
because there are atmospheric layers moving in opposite di- 
rection [i.e. the atmosphere is not in free fall]. The approxi- 
mation (|6]) is not valid in these phases in spite of a' ang ' ~ 
for a short time when the rapidly changing a^ ang - ) (r, t) has 
a sign change. These shocked phases were therefore not in- 
cluded in i = 1, 2, n. 

The right hand side of Eq. consists of quantities de- 
rived directly from the photometry without differentiations 
of the polynomial fits. This can be determined with an ac- 
curacy of < 0.15 from the UBV(RI)c photometry in each 
k. 

(iii) By introducing the averaged value 



■Ma 

d 2 



n ^ 

i=i 



d 2 



G 



g c (t) - a(r,t) 



Jdyn) 



(r,t) 



(4) 



we reduce Eq. ((4]) to an algebraic equation for unknown dj 
at any tj, j = 1,...,N if a' dyn ' = 0; the largest term is 
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linear in dj . These -/V equations can be solved by elementary 
operations for all j. The distance was obtained from 



iV 



(8) 



At this step, an averaged dynamical correction 



a (^)(R,N) 



3 = 1 

JV 



(dyn) (7?,^) 



j'=i 



a(R, tj 



(9) 



was defined with M a d~ 2 and d(iV) from Eqs. Q and ((Sj 
facilitating the reduction of the number of the equations Q 
to be solved. 

(iv) To obtain the final d(N (u ^) the step (iii) was 
repeated, reducing N. One or two iterative steps 
were necessary to obtain a( dyn )(ii, iV*- 11 ') <C 



and 



,(dyn) 



.[«(<)], 

=!,...< 0.3flb[«(*)] where ffl.[fl(t)] = 



GMaR~ 2 (t). That is, the roots dj of Eq. (0}, originating 
from the tj intervals, were removed from Eq. (J5J) which pro- 
duced outlier values with a (dyn) (_R, tj) > 0.3g s [R(tj)]. The re- 
sults of the iterative steps are given in Table [4] 

The interval of the search for roots of Eq. Q was limited 



to 



(10) 



To obtain this upper limit, we took into account only the 
linear term of the acceleration when the atmosphere is hit 
by the strongest shock wave; this is the state of maximal 
compression. At this epoch the atmosphere is just reversed 
from inward to outward motion. That is, v(dv/dr) ~ 0, 
dv/dt — #d > 0, R is minimal, the deceleration by the 
static gravity is negligible (g s <C dv/dt) and, furthermore, 
a (dyn) ^ q wag SU pp 0se( j However, the upper limit d m ax was 
taken into account only if the satisfaction of these conditions 
plus, C (I) and C (II) , was verified afterwards. 



the different colour-colour combinations and the scatter of 
i? when it is determined from the brightness in V, Rc and 
V plus bolometric correction and the Stefan-Boltzmann 
law. The components of the error of L cq are composed of 
errors from (# 2 T 4 ), M^dT 2 and d, respectively. The nu- 
merical values are ±0.18, ±0.35, ±0.70 for GSC 4868-0831 
±0.2, ±1.3, ±3.7 for V372 Ser. A summary of the parame- 
ters is given in Table [5] 

The error of the magnitude averaged (My) was esti- 
mated from sum of the errors of M a d~ 2 and d in Table [J 

Fig. presents an insight into the variable quantities 
of the atmosphere at a shock-free and a shocked phase for 
both stars. The segment of V(t) is plotted in the uppermost 
row for orientation. The scatter of T e (t) (row 2) reflects the 
scatter in our photometry. R(t) and the velocities and accel- 
erations were computed from the smoothed $(t) and log g c (t) 
using d and A1 a from Table [4] 

We see a considerable gradient of v(r) for GSC 4868- 
0831 at HJD > 2454829.6, just when a descending branch 
of the light curve ended. The deeper layers were moving 
outwards much faster because the atmosphere was hit by a 
shock wave, and the velocity excess reached as 25 kms . 
Otherwise \v(R) — v(R — hg 1 )\ ~ a few km s _1 was found 
for both stars. 

Row 5 is a plot of g s (R,t). The different types of 
GSC 4868-0831 and V372 Ser are obvious if we compare 
R(t) and g B (R,t); GSC 4868-0831 is a subgiant star while 
V372 Ser has a true RRd character. The components of 
the acceleration are plotted in row 6 in the relative units 



<7d = a 



(dyn) 



(R,t)g~ 1 (R,t), q r = R(R,t)g7\R,t), q a = 
a{R,t)gJ 1 {R,t). Here, |qd | < 0.3 is the domain when C"' 
is satisfied if the light curve is not in ascending branch. 
The rapid changes of the components were not smoothed 
out, especially those of q d of V372 Ser at HJD> 2454250.57. 
These originate partly from violating C' 1 ', and partly from 
the stronger effect of the thermic shock on the more dilute 
atmosphere of V372 Ser. qd > 10 means that g c derived from 
the photometry has a very loose connection with the atmo- 
spheric kinematics. This must be corrected by a term 3> g s 
to obtain the actual acceleration in the atmosphere. 



4.1 Numerical results 

Our observations provided 10 light curve segments of suf- 
ficient length and quality to calculate the polynomial fits 
of log g e (t) , T e (t) and i?(£) for both stars. The order of 
magnitude limits were found to be — 10 -20 < 0[a' ang '] < 
10 _18 rad-s -2 for both stars. The acceleration-free intervals 
can be found for Eq. Q from quantities that are of the form 

/W±o(d- 1 ). 

The results from step (i), [M] and E(B — V) are given 
in Table [3] The results from steps (ii)-(iv) and the final d 
and Ma. are summarized in Table [4] The asterisk in the last 
column denotes the value of iV' 11 '. 

The position of the stars in a theoretical HRD was cal- 
culated by 

T cq = (i9 2 r e 4 ) 1 ' /4 (i9) _1//2 , (11) 
L cq = 4W 2 <tf 2 (v9)T < 4 ( V 9)) (12) 

l|Carnev. Strom fc Jones! I1992D . The numerical values are 
given in Table [5] The remarkably small error of T eq orig- 
inates from the sum of the scatter of T , determined from 



5 DISCUSSION 

5.1 Remarks on the photometric inverse problem 

Because of the lack of high-dispersion spectra the relations 

P = /(CI 1 ,CI 2 ...) (13) 

play a crucial role. Here, V and Cli denote an astrophysical 
quantity and colour index, respectively. For example, a par- 
ticular form of / belongs to V — T c , V — log g , etc. These ar e 
given in tabular form for the ATLAS models (|Kuruczlll997h , 
and their parameters are E(B — V), [M], etc. (As emphasized 
in Sec. 3, [M] is on the solar scale and / would have another 
form for a peculiar atmospheric chemical composition.) 

The derived fundamental parameters depend on the 
reddening and metallicity. A larger E(B — V) leads to larger 
log g and T c , the dependence on [M] is similar, but to a lesser 
extent. The essence of our method is that we search for the 
minimal standard errors A log g c , AT C of log g c , T c from some 
30 colour index pairs as a function of [M] and E(B — V). 
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Table 4. The results from steps (i)-(iv). 



{E{B — V), [M]} 


n 


M & d~ 2 x 10 7 d 




a (dyn) ( Ri t ) 


N 


[mag] , [dex] 






[pc] 


[M Q ] 


[ms -2 ] 










GSC 4868-0831 








{0.008, -1.05} 


7 


40.3 ±6.7 


451 ± 99 


.82 ± .36 


-58.4 


280 








467 ± 16 


.88 ± .06 


0.10 


21* 








V372 Ser 








{0.003, -0.53} 


20 


6.12 ± .31 


1183 ± 74 


0.85 ± .13 


13.90 


529 








1020 ± 87 


0.63 ± .12 


1.64 


249 








964 ± 81 


0.57 ± .10 


0.41 


135* 



Note. The * indicates N^; |a( d y n >| < 3ms" 2 is satisfied for all tj,j = 1, JV' 11 ). 




Figure 3. The variable physical parameters of GSC 4868-0831 and V372 Ser in descending and ascending branches of the V light curve. 
Row 1: light curves V(t). Row 2: T c (t). Row 3: R(t). Row 4: solid lines: v(R,i), dotted lines: v(R - h^ 1 ,^. Row 5: g s (R,t). Row 6: 
relative accelerations (solid lines: = a^- dyn \R,t)g7f 1 (R,t), dotted lines: <jy = R(R, tygTf 1 (R,t), dashed lines: q a = a(R,t)gJ 1 (R,t)). 



This method is self-consistent to the highest degree from an 
astrophysical point of view. We use all colour information of 
the atmospheric models to determine the relevant parame- 
ters and quantities. In doing so, we can avoid the systematic 
errors inherent to an arbitrary choice of relatio ns such as, for 
exam ple, T e -(V - I c ), [M/H]-(B - V), (e.g. iDekanv et all 
2008). In addition, we do not use semi-empirical calibra- 
tions at all (e.g. the T -colour index relations derived for 
non-variable main sequence stars of normal chemical com- 



position and, as next step, g eneralized to variable g iant stars 
with large metal deficiency. IClementini et alll2000l ). 

Of course, in comparison with a line by line analysis 
of high-dispersion spectroscopy, a multicolour photometry 
is sensitive only to the general shape of the optical flux 
of the star. However, it is possible to find the model at- 
mosphere by reproducing the continuum flux in the visi- 
ble wavelength interval. Our method can provide a global 
parameter of the chemical composition in the atmosphere 
summarized as metallicity [M]. This [M] is the most suit- 
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able for our purposes, because it accounts for the effect of 
all elements on the optical continuum, including those that 
do not have lines. Of course, it might differ from the overall 
metallicity of the star (used in the theory of stellar struc- 
ture, as the pulsation does not stir up the deepest layers 
where the nuclear reactions take place) or from the aver- 
aged metallicity derived from high dispersion spectroscopy. 

The fields are at high Galactic latitude: GSC 4868- 
0831, b = +24°; V372 Ser, b = +45°. Their redden- 
ing is very small, and therefore the parameters (logi? c , 
T c , d, Ma) reported in Section 4 are lower limits at the 
same time if [M] is fixed. The upper limits are E(B — 
V) < 0.023, 0.085 dSchlegel. Finkbeiner. fc Daviej|l998T ). or 
E(B-V) < 0.03.0 045 lBurstein fc Heiledll982l , "respectively. 
We attribute our smaller E(B — V) to two factors: our 
method measures the reddening of a point source and the 
excess reddening must originate from a region beyond GSC 
4868-083f , V372 Ser, and the comparison stars. 

The reddening E(B - V) = 0.085 of V372 Ser derived 
from the D1RBE differs from our value above the 5a level. 
Taking this large E(B — V) would result in an unacceptable 
increment from (AT e ) 2 49 = 6.2 + 0.2 to 8.2 + 0.4. The param- 
eters would increase to logg c = 4.11, T e = 7092 leading to a 
mass above 5Mq which cannot be reconciled with any ac- 
tual theoretical knowledge about pulsating stars. Therefore, 
E(B -V)> 0.006 can be ruled out. 



log L a 




Figure 4. The positions of the stars in a theoretical HRD. The 
square denotes SU Dra, (with [L e q,T eq ] were taken from Paper 
I), The triangle and circle denote V372 Ser and GSC 4868-0831, 
respectively. The solid lines from left to right denote the blue and 
red edges of the instability str ip. The dashed line denotes the zero 
age horizontal branch of M3 iSilva Aguierre et al.ll200gh . 



5.2 The use of U observations 

Concerning GSC4868-0831, an interesting result can be seen 
from Fig. [2] the phases satisfying and not satisfying C^ 1 ' seg- 
regate clearly in the colour-colour diagrams and completely 
different T e (t) and log g e (t) are obtained for these phases if a 
BVIc photometry only is used as input. Fig. [^demonstrates 
that reliable T e (t) and g e {i) can be determined only if colour 
indices containing U are used in the T e , logg domain of RR 
stars. The use of one colour index or BVIc photometry only 
is not sufficient and can be misleading, even if it is limited 
to determining T c only. 

Furthermore, our U observations might reveal the sub- 
giant charac ter of GSC 4868-083 1 which was not suspected 
previously dWils fc Oterol 120051. I Wils. Lloyd fc Bernhard] 



2006, ISzczygiel fe Fabrvczkv 20071 ') . An important conclu 
sion has emerged that it is not possible to determine the 
luminosity class of a pulsating star if only periods or period 
ratio are available. Multicolour observations, covering the 
ultraviolet, are needed to classify DM pulsators properly. 

5.3 Remarks on the QSAA 

We have to emphasize that L eq and T eq are first approxima- 
tions from Eqs. (|11[) and ()12[) because some elements of the 
averaging were obtained assuming QSAA in phases when 
C (D 

was violated. Qualitative considerations suggest a posi- 
tive correction to T e (t) of QSAA in the shocked phases when 
excess radiation and dissipation exist from shock waves. 
Corrections emerging from a dynamical model atmosphere 
would not modify the main fundamental parameters M a 
and d because they were determined from phases when both 
quantitative conditions of the validity of QSAA were satis- 
fied. Ma, and d can be considered as well substantiated em- 
pirical data from the ATLAS static model atmospheres plus 



some basic hydrodynamics. Dynamical model atmospheres 
are beyond the scope of this series of papers. 

The sampling of the quasi-repetitive curves introduced 
negligible error, which can be esti mated by comparing T e 
and i? from the fitted colour curves (|Benko fc Barczall2009h 
with those from the N = 529 observations of V372 Ser. 

We remark that d = (1145 ± 73) pc, Ma, = (0.83 ± 
.17)M e are the results for {E(B - V) = 0.003, [M] = 
—0.53} if the UAA is applied, that is, if dv/dr = is as- 
sumed in Eq. Q (and, as a consequence, i9 = R/d and 
dv/dt — dd). This distance and the change of _M a d~ 2 = 



6.12 ± .31 ->■ 7.64 ± .82 gives L, 



(UAA) 



= 42.8Z/Q which put 



V372 Ser i n a position just at the lower lim it of stable DM 
pulsation (|Szab6. Kollath fc Buchleil 12004). Of course, the 
physical input of the UAA is much less than that of our 
extended hydrodynamic treatment represented by Eqs. (|4]), 
{5). fn spite of the better agreement, the data from UAA 
must not be accepted because the UAA is a rigid and less 
realistic approximation in comparison with a compressible 
model atmosphere. 



5.4 Classification of the stars 



The positions of our stars and SU Dra in a theoretical HRD, 
(i.e. [L eq , T eq ]) are plotted in Fig. U For orientation, the 
instability strip and the zero age horizontal branch of the 
metal defici ent ([M]« -1.6) globular cl uster M3 (NGC 5272) 
are shown l|Silva Aguierre et al.ll2008T ). We emphasize that 
[I/ e q,7eq] correspond to L and T e of non- variable stars, and 
they can be directly compared with those from the theo- 
retical studies on stellar structure, pulsation and evolution. 
The only source of error is the difference of T e (t) and i9(t) 
in Eqs. (|lll 11 2 p from QSAA and dynamical model atmo- 



Fundamental parameters of RR Lyrae stars. II 9 



spheres, respectively. An error has not been propagated into 
the position of the star by semi-empirical relations like T c - 
(B - V) etc. 

The period ratios Pi/Pq are = 0.7443 and 0.7450 for 
V372 Ser and GSC 4868-0831, respectively. These are in the 
canonical range of RRd stars, as Ma of V372 Ser is also. 
The equilibrium effective temperatures of both stars are in 
the 5630 < T cq < 7080 K interval where DM pulsation can 
be stable. The results of our analysis confirm that V372 Ser 
is a true RRd star. However, it is subluminous by a factor 
of « 2-3 in comparison with the puls ation models of RRd 
stars l|Szabo, Kollath fe Buchleill2004r ). 

Ma of GSC 4868-0831 exceeds somewhat the canonical 
mass of RRd stars; its subluminosity is even larger, ~ 4-5 
and its radius is less than half that of an RR star. The 
periods exclude its identification as DM SX Phe star. To 
the best of our knowledge, GSC 4868-0831 is the first known 
subgiant star pulsating in two modes of large amplitude. It 
is slightly (~ 1 mag) above the zero age main sequence of 
normal chemical composition, and it is below the extension 
of the instability strip of M3. 



6 CONCLUSIONS 

We have reported 302 UBV{RI)c observations of GSC 
4868-0831. To the best of our knowledge, this is the sec- 
ond DM pulsating star with well documented photometric 
behaviour in the optical and near ultraviolet bands. 

We have determined, for the first time, fundamental pa- 
rameters of DM pulsators using only the theory of stellar at- 
mospheres. Our method is purely photometric; the eventual 
uncertainty of observing radial velocities and their conver- 
sion into the reference frame in the centre of the star can be 
avoided using our method. Our method has made extensive 
use the ATLAS model atmospheres and some basic hydrody- 
namics. Our astrophysical input is completely different from 
the theory of stellar pulsation and evolution. Therefore, our 
method can be used as a check or challenge for these more 
involved theories. 

We have summarized our results in Table [5] We have 
found the surprising result that GSC 4868-0831 is not an 
RRd star, but a subgiant star pulsating in two modes. We 
have found a significant subluminosity of the RRd star V372 
Ser in comparison with the luminosity from present day the- 
ory of stellar pulsation. These objects must be considered 
as a challenge to extend the search for stable DM pulsation. 
Qualitative considerations suggest that dynamical model at- 
mospheres would have higher T oq and L cq ; both corrections 
would shift the points in Fig. [4] upwards and to the left with 
respect to their position from QSAA. However, we must 
not forget that our results from static ATLAS model at- 
mospheres have yielded the ratio Ma,d~ 2 and d from phases 
when a dynamical model atmosphere is not necessary (i.e. 
QSAA is a reliable approximation of the pulsating atmo- 
sphere). We are faced with the dilemma that there is either 
subluminosity with acceptable mass or acceptable luminos- 
ity with too large mass if we are to reconcile these parame- 
ters from the present-day theory of stellar pulsation. 



Table 5. Summary of the parameters. 



GSC 4868-0831 V372 Scr 



Po 


0^5649* 


ff?4712289* 


A 


0.07271* mag 


0.1534* mag 


Pi 


0? 42085* 


0'?3507310t 


Ai 


0.15815* mag 


0.2059* mag 


E(B - V) 


0.008 ± .002 mag 


0.003 ± .003 mag 


[M] 


-1.05 ± .10 dex 


-0.53 ± .05 dex 


d 


467 ± 16 pc 


964 ± 81 pc 


(M v ) 


+2.40 ± .24 mag 


+1.58 ± .22 mag 


-Ma 


0.88 ± .OGMq 


0.57 ± .10M & 


-^min 


1.97R 


4.O7iJ 




2.05i? Q 


4.4Oi? 


Tcq 


6924 ± 35K 


6722 ± 20K 


^cq 


8.97 ± 1.23L 


21.9 ± 5.2L 



Notes. 

*: From lBenko fe Barczal [|2009ri . 

*: determined from the V observations. 
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APPENDIX A: DESCRIPTION OF THE 
PROGRAM PACKAGE BBK 

The package BBK is designed to extract the fundamental 
parameters of RR stars from high-quality UBV(RI)c ob- 
servations. A brief description and some technical details 
are given here. 

It is important to have the colours and colour indices 
as close as poss ible to the inte rnational colour system which 
was applied bv lKurucd (|1997T ) in order to convert the phys- 
ical fluxes of the ATLAS models to the stellar magnitude 
and colour system. An error (or errors) in the colour indices 
can lead to false results. To obtain reliable results, some 
N > 300 five colour observations are needed. These must be 
distributed uniformly over a representative light curve be- 
cause the physical quantities if, ho must be differentiated. 

Proper use of the program package requires basic knowl- 
edge of absolute stellar photometry, the theory of stellar 
atmospheres and hydrodynamics. As iKurucj (1997) said: 
'Neither the programs nor data are black boxes. You should 
not be using them if you do not have some understanding 
of the physics and of the programming in the source code.' 
This warning is appropriate for BBK. 

BBK consists of FORTRAN p rograms, z-sh ell scripts, 
input files, data files, and a manual (Barcza 2011). UNIX or 
LINUX environment, installation of a FORTRAN compiler, 
z-shell, graphical packages (e.g. SUPERMONGO, or GNU- 
PLOT) are necessary. Five steps are involved, a detailed de- 
scription of which can be found in the manual. Inspection, 
evaluation, plots of the partial results are necessary before 
continuing to the next step. 

Step (I) The atmospheric metallicity and the interstellar red- 
dening of the star are determined from selected phases when 
the atmosphere is free of shocks (i.e. from colour indices in 
the descending branch). 

Step (II) The conversions 



/(CIi,CI 2 ...) 



(Al) 



are executed to select the static ATLAS models with the 
best fit to the observed and theoretical colour indices to 
obtain log g e and T e as a function of phase. The variation of 
•& is determined by comparing the physical fluxes of the star 
with those of the selected theoretical models. The physical 
fluxes V, Rc, V+ bolometric correction are used, those of 
the star are calculated from the absolute calibration of Vega 
(|Tiig. White fc Lockwoodlll977i U 

Step(III) Polynomial fits are calculated to obtain log<? e , T e , 
ho, angular velocity and acceleration as a function of 



phase. The upper limit dj < d n 



must be fixed 



by inspecting ■& as a function of phase. 

Step(IV) The fits are introduced in the Euler equation of 
hydrodynamics to determine the angular acceleration and 
to find the acceleration-free phases of the atmosphere. The 



transient acceleration-free intervals in the ascending branch 
(i.e. in the shocked phases) must be manually deleted from 
the averaging to find M^d~ 2 . 

Step(V) •&, ho, dho/dt, g e and Ma,d~ 2 are introduced 

in Eqs. (I4I5P and Eq. @ is solved for d with the whole (or 
partial) set containing N (or < N) phase points. The upper 
limit a (dyn ' uppcr limit) must be specified to exclude the phase 
points from the final solution for d [i.e. from A^ 11 -*] to achieve 
in all phase points giving d. One 



,(dyn)| 



< a 



(dyn, upper limit) 



or two runs might be necessary with decreasing TV. 



